Exploring the spin-| frustrated square lattice model 
with high-field magnetization measurements 
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We report on high-field magnetization measurements for a number of layered vanadium phosphates 
that were recently recognized as spin- 1/2 frustrated square lattice compounds with ferromagnetic 
nearest-neighbor couplings (Ji) and antiferromagnetic next-nearest-neighbor couplings (J2). The 
saturation fields of the materials lie in the range from 4 to 24 T and show excellent agreement with the 
previous estimates of the exchange couplings deduced from low-field thermodynamic measurements. 
The consistency of the high-field data with the regular frustrated square lattice model provides 
experimental evidence for a weak impact of spatial anisotropy on the nearest-neighbor couplings 
in layered vanadium phosphates. The variation of the J2/J1 ratio within the compound family 
facilitates the experimental access to the evolution of the magnetization curve upon the change of 
the frustration magnitude. Our results support the recent theoretical prediction by Thalmeier et al. 
[Phys. Rev. B, 77, 104441 (2008)] and give evidence for the enhanced bending of the magnetization 
curves due to the increasing frustration of the underlying spin system. 



PACS numbers: 75.50.-y, 75.10.Jm, 75.30.Et 

Quantum phenomena and exotic ground states in low- 
dimensional magnets receive broad attention from both 
experimental and theoretical sides. The application of 
an external magnetic field provides a natural and ex- 
perimentally feasible way to alter the properties of such 
systems. Recent studies showed that numerous low- 
dimensional magnets undergo field-induced phase tran- 
sitions that can be described with bosonic models^ A 
case of particular interest are the frustrated spin sys- 
tems, because frustration often leads to the localization of 
magnetic excitations on the lattice and can cause the for- 
mation of unusual ordered phases, evidenced by plateau 
features on the magnetization curves. 

The frustrated square lattice (FSL) is one of the most 
simple low-dimensional frustrated spin models. This 
model assumes the competing interactions along the sides 
and the diagonals of the square, J\ and J2, respec- 
tively (see Fig. [T]). The variation of the frustration ra- 
tio a = J2/J1 enables to tune the system from a non- 
frustrated to a strongly frustrated regime. For the case of 
spin-i, theoretical studies suggest three ordered phases 
and two critical spin-liquid regions in the ground state 
phase diagram.-^* 3 . Model simulations propose the forma- 
tion of a magnetization plateau at half-saturation^ (and, 
possibly, another plateau at one-third of the saturated 
magnetization^) for the region of J2/J1 — 0.5. Later 
theoretical studies focused on the magnetocaloric effect, 
the saturation fields, and the shape of the magnetiza- 
tion curves for the whole phase diagram. 6 ' 7 In particu- 
lar, frustration-driven quantum fluctuations were shown 
to enhance the bending of the magnetization curves as 
compared to the linear curve for the classical system* 7 - 
Despite the broad theoretical interest, systematic exper- 
imental high-field studies of FSL compounds have not 
been reported. 



Recently, a number of spin-i FSL model compounds 
were proposed i 8 ' 9 ' 10 ! 11 ! 12 Most of these materials reveal 
relatively weak exchange couplings, thus efficient high- 
field experiments in pulsed or even static fields are fea- 
sible. In this Brief Report, we will focus on the family 
of layered vanadium phosphates that lie in the region of 
columnar antiferromagnetic (CAF) ordering in the ferro- 
magnetic (FM) Ji - antiferromagnetic (AFM) J2 part of 
the general phase diagram (see Fig. Q]). Below, we will 
discuss five compounds of this family. The chemical com- 
positions of these compounds are listed in Table Q] along 
with their abbreviations. 

The low-field properties of PVO and SZVO have been 
reported in Refs. H and[H. The compounds NVOF and 
PZVO were recognized as the FSL materials quite re- 
cently, and the detailed study of their low-field behavior 
will be published elsewhere. The investigation of BCVO 
has been reported in Ref. along with the full magneti- 
zation curve that will be used here for the sake of compar- 
ison. The high-field data for the other four compounds 
are original and have not been previously reported. 

Layered vanadium phosphates have a common struc- 
tural feature, the [VOPO4] magnetic layers with V +4 C>5 
square pyramids and superexchange pathways for the 
nearest-neighbor (NN) and next-nearest-neighbor (NNN) 
interactions (see Fig. [I]). The low crystal symmetry usu- 
ally leads to the non-equivalent couplings along the sides 
(Ji, J[) and the diagonals (J2, J 2) of the square, thus a 
rigorous spin model including first and second neighbors 
should comprise four independent parameters. However, 
our recent computational studyi^ demonstrated that this 
effect of the spatial anisotropy is usually weak. In partic- 
ular, thermodynamic properties of the layered vanadium 
phosphates can be accurately described with the regular 
FSL model and do not allow to estimate the magnitude 
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TABLE I: Averaged NN (Jj) and NNN (J 2 ) exchange couplings in layered vanadium phosphates according to the low-field 
thermodynamic measurements, the frustration ratio a, thermodynamic energy scale J c = y Jf + Jf , the experimental satura- 
tion field (noH s ), and the calculated saturation field for the regular FSL model (/xo-ffJ eg ), as determined from Eq. ^ using the 
Ji and J 2 values. 



Index 


Compounds 


Jx (K) 


h (K) 


a = J 2 /Ji 


Jc (K) 


» H S (T) 


Mo#r s (t) 


NVOF 


Na 1 . 5 VOP0 4 Fo. 5 


-3.2 


6.3 


-2.0 


7.1 


15.4 


14.3 


PZVO 


PbZnVO(P0 4 ) 2 


-5.2 


10.0 


-1.9 


11.3 


23.4 


22.6 


PVO 


Pb 2 VO(P0 4 ) 2 


-5.1 


9.4 


-1.8 


10.7 


20.9 


20.9 


SZVO 


SrZnVO(P0 4 ) 2 


-8.3 


8.9 


-1.1 


12.2 


14.2 


14.5 


BCVO (Ref. 14) 


BaCdVO(PQ 4 ) 2 


-3.6 


3.2 


-0.9 


4.8 


4.2 


4.3 




FIG. 1: Left panel: part of the ground state phase diagram 
for the FSL model^ solid filling shows the regions of the FM 
and CAF ordering, while hatched filling denotes the critical 
region around a — —0.5. Right panel: the [VOP0 4 ] layer in 
the structures of layered vanadium phosphates. In the left 
panel, the model compounds under discussion are placed into 
the phase diagram according to Table [I] 



of the spatial anisotropy. In contrast, the field depen- 
dence of the magnetization should be altered by the spa- 
tial anisotropy of nearest-neighbor couplings. This prob- 
lem will be further discussed below. 

The low-field thermodynamic data yield averaged 
nearest-neighbor and next-nearest-neighbor couplings 
[Ji = (Ji + J0/2 and J 2 = (J2 + J^/2, respectively] for 
the distorted square lattice.— The resulting estimates of 
Ji and J2 are listed in Table fl] These numbers can be 
used to calculate the effective frustration ratio a — J 2 / Ji 
and to place the compounds into the general phase dia- 
gram (Fig. [[]) . 

Polycrystalline samples of PVO and SZVO were pre- 
pared according to Ref. 1 1 31 The powder samples of 
NVOF and PZVO were obtained by firing the mixtures 
of NaP0 3 , NaF, and V0 2 (NVOF) or PbZnP 2 7 , V 2 3 , 
and V 2 5 (PZVO) in evacuated silica tubes at 700 °C. 
The details of the preparation procedure will be reported 
elsewhere. The phase purity was checked using X-ray 
diffraction. All the samples except for PZVO were single- 
phase. The sample of PZVO contained 2 — 3% of the 
diamagnetic PbZnP 2 07 impurity that did not influence 
on the shape of the magnetization curve. 

Magnetization measurements were performed in the 
Dresden High Magnetic Field Laboratory using a pulsed 



magnet powered by an 1.44 MJ capacitor bank. With 
an inner bore of 20 mm, the magnet yielded fields up to 
60 Tesla with a rise time of 7 ms and the total pulse du- 
ration of about 20 ms. Magnetic moment of the sample 
was obtained by integration of the voltage induced in a 
compensated pick-up coil system surrounding the sam- 
ple. Each measurement was performed twice, with and 
without the sample, in order to remove the background 
signal contribution. All the data were collected at a tem- 
perature of 1.4 K. The curves measured on increasing and 
decreasing field coincided, thus indicating the lack of any 
irreversible effects upon magnetization and demagnetiza- 
tion of the samples. 

All the four compounds show monotonous increase of 
the magnetization (M) upon the application of the ex- 
ternal magnetic field (H). The saturation is evidenced 
by a sharp bend of the curve (see Fig. [2]). Above this 
anomaly, the magnetization reaches the saturated value 
(M s ) and remains field- independent. 

The saturation field (H s ) is usually determined by 
taking derivatives of finite-temperature magnetization 
curves. The H s value is assigned to the inflection point 
of the first derivative or to the minimum of the second 
derivative (both refer to the point of the maximum curva- 
ture for the M vs. H dependence) . The saturation fields 
obtained in this way are listed in Tabled] The uncertainty 
of the resulting numbers depends on two factors: (i) the 
accuracy of the experimental data; (ii) the temperature 
of the measurement. Experiments in pulsed fields yield a 
very high field resolution; therefore, the calculation of the 
derivatives is quite accurate, and the uncertainty from (i) 
does not exceed 0.1 T. The effect of (ii) is a smoothen- 
ing of the saturation anomaly. At finite temperatures, 
the maximum magnetization is reached above H s , be- 
cause an additional energy (and, hence, higher field) is 
required to suppress thermal fluctuations. 

To illustrate the finite-temperature effect and to es- 
timate the resulting error in the H s values, we sim- 
ulated magnetization curves for a representative two- 
dimensional spin model, the unfrustrated square lattice 
with the single parameter J — J\ (J2 = 0). The simu- 
lations were performed using the ALPS code 16 and the 
quantum Monte-Carlo (QMC) directed loop algorithm 
in the stochastic series expansion representation. 17 The 
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FIG. 2: Experimental magnetization curves normalized to 
the saturation field {H s ) and to the saturated magnetization 
(M s ), the arrow shows the evolution of the curves upon the 
increase of the frustration ratio a = J2/J1. The inset shows 
QMC simulations for the unfrustrated square lattice with the 
exchange coupling J. Solid, dashed, and dash-dotted curves 
refer to the temperatures (T / J) of 0.02, 0.1, and 0.2. Solid ar- 
rows denote the points of the maximum curvature, while the 
dashed arrow shows the shift of these points with the increase 
of temperature. 



QMC magnetization curves for the unfrustrated square 
lattice are shown in the inset of Fig. O One can see that 
thermal fluctuations tend to smoothen the curve, and 
the saturated value M s is reached well above H s . Still, 
the points of the maximum curvature remain close to the 
true saturation field for T = 0, although the experimental 
H s numbers at finite T should be slightly overestimated. 
For example, the error amounts to about 1% of H s at 
T = 0.1J and to about 2% of H s at T = 0.2J. For the 
frustrated square lattice, the energy scale is defined by 
the value of J c = \J J 2 + Jf . The experiments for PZVO, 
PVO, and SZVO are performed at T ~ 0.1J C , thus the 
uncertainty of H s amounts to 0.15 — 0.20 T. In case of 
NVOF, the uncertainty may be slightly higher due to the 
lower J c . The magnetization curve of BCVO has been 
collected at 0.5 K— that also corresponds to T ~ 0.1J C . 
Overall, the total error from (i) and (ii) should not exceed 
0.3 T for all the compounds. 

After establishing the uncertainty of the experimental 
H s values, one can compare them to the low-field esti- 
mates of the exchange couplings. All the compounds un- 
der investigation lie in the region of CAF ordering (see 
Fig. [TJ, hence the saturation field of the regular FSL 
amounts tc 6 



fi H s = 45(Ji + 2J 2 )k B /{gHB), 



(1) 



where S = \ is the spin value, k B is the Boltzmann 
constant, fi B is Bohr magneton, and g is the (/-factor. 
In the following, we assume g = 1.95 as a representative 
value for V +4 -containing compounds (see, e.g., Ref. [l8[) . 



The experimental H s values and the low-field estimates 
are compared in Table[I] For most of the compounds, we 
find perfect agreement with the previously obtained Ji 
and J2 values. In case of NVOF and PZVO, the exchange 
integrals slightly underestimate H s , and this effect will 
be further discussed below. A comparison of the satu- 
ration fields manifests the role of the frustration in lay- 
ered vanadium phosphates. The saturation field of SZVO 
is well below that for PZVO and PVO, although SZVO 
shows the largest J c among the whole family. This con- 
clusion is consistent with the previous phenomenological 
evidence of the frustration, e.g., the suppression of the 
specific heat maximum in SZVO as compared to PVO^ 

The saturation fields can also be used to estimate the 
spatial anisotropy of nearest-neighbor couplings in lay- 
ered vanadium phosphates. Above, we have mentioned 
that the low crystal symmetry leads to incquivalent ex- 
change couplings along the sides and the diagonals of the 
square (Fig. [l}. In case of the CAF ordering, the differ- 
ence between J\ and J[ should alter the saturation field. 
If J 1 and J[ are FM, and | J[ \ < \Ji\, the columnar order- 
ing is favored by J\ , J2 , and J' 2 , while disfavored by J[ . 
Thus, the reduction of |J{| should stabilize the columnar 
ordering with respect to the ferromagnetic state, thus 
leading to the increase of H s . Using the classical Heisen- 
berg model, one easily finds that the increase of the true 
saturation field (H s ) as compared to the estimate for 
the regular lattice (Hl cs ) is proportional to (J[ — Ji): 
fi Q H s - (i H™z = 2S(J[ - Ji)k B /(gfi B ). For PVO, 
SZVO, and BCVO, the H s and iJ s reg values nearly match. 
Yet PZVO shows the anisotropy J[ - J x ~ 0.8 K. In 
case of NVOF, the anisotropy is even more pronounced: 
,/{ — Ji ~ 1.4 K. The origin of the spatial anisotropy in 
these compounds is still to be clarified. 

Now, we turn to the shapes of the experimental magne- 
tization curves. Fig. [2] presents a comparative plot with 
all the curves normalized to the saturation field and to 
the saturated magnetization. To determine M s , we used 
M values well above H s . One can observe a clear differ- 
ence between the curves. As the frustration is increased, 
the bending gets more pronounced. The strongest bend- 
ing is found for BCVO with the largest frustration ratio 
a and the strongest frustration. This result confirms the 
recent theoretical prediction on the quantum effects in 
the high-field magnetization of the FSL models Quali- 
tatively, the effect can be understood as the competition 
of quantum fluctuations and the external field. Frus- 
tration induces strong quantum fluctuations that reduce 
magnetization values at low fields and cause the positive 
curvature for the field dependence of the magnetization. 

The agreement between experiment and theory is not 
only qualitative, but quantitative as well. In Fig. [31 
we plot the experimental data for PZVO along with 
the simulated magnetization curve. The simulations 
were done for clusters with 16, 20, and 24 sites using 
finite-temperature Lanczos method (FTLM)^ The data 
points were obtained via the Bonner-Fisher technique. 20 
The frustration ratio was taken from the experimental es- 
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FIG. 3: Experimental magnetization curve for 
PbZnVO(P0 4 )2 along with the FTLM simulations. Circles, 
squares, and diamonds denote the results for the 16-site, 
20-site, and 24-site clusters, respectively. 

timate J2/J1 = —1.9. We find perfect agreement between 
the experimental and simulated curves above 0AH s . At 
lower fields, the simulated curve lies below the experi- 
mental one. This should be related to finite-size effects 
which are expected to be most pronounced at low fields. 
Nevertheless, the consistency of the experimental and 
theoretical results shows that simulations for the rela- 
tively small clusters can be an effective tool for studying 
high-held magnetization of different frustrated spin sys- 
tems in two dimensions (e.g., the kagome lattice). 

In conclusion, we performed high-field magnetization 



measurements for a number of layered vanadium phos- 
phates representing spin-i FSL compounds. The pulsed- 
field experiments provide accurate estimates of the satu- 
ration holds, while the latter hold important quantitative 
information on the properties of the FSL systems. The 
saturation fields of layered vanadium phosphates are con- 
sistent with the previous low-field estimates of averaged 
exchange couplings J\ and J2. This confirms the place- 
ment of the compounds on the phase diagram and evi- 
dences the strongest frustration in SrZnVO(P04)2 and 
BaCdVO(P04)2- Additionally, the values of the satura- 
tion fields provide an estimate for the spatial anisotropy 
of nearest-neighbor couplings. We show that this effect 
is relatively weak, and the difference between J[ and Ji 
does not exceed 1.0 — 1.5 K in the whole family of layered 
vanadium phosphates. The comparison of the magnetiza- 
tion curves for the compounds with different frustration 
ratios reveals the effect of frustration-induced quantum 
fluctuations. The frustration leads to the increased bend- 
ing of the magnetization curves in agreement with the 
theoretical predictions. The direct comparison of the ex- 
perimental and simulated curves indicates high accuracy 
of the simulation techniques and validates further model 
studies of two-dimensional frustrated spin systems. 
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